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Abstract 

We apply robust control technics to an adaptive optics system including a dynamic 
model of the deformable mirror. The dynamic model of the mirror is a modification 
of the usual plate equation. We propose also a state-space approach to model the tur- 
bulent phase. A continuous time control of our model is suggested taking into account 
the frequential behavior of the turbulent phase. An Hoa controller is designed in an 
infinite dimensional setting. Due to the multivariable nature of the control problem 
involved in adaptive optics systems, a significant improvement is obtained with respect 
to traditional single input single output methods. 

Keywords: Adaptive Optics, Robust control, Partial differential equations. 
OCIS: 010.1330, 120.4640. 



1 Introduction 

For several decades it has been now possible to use adaptive optic (AO) systems to ac- 
tively correct the distortions affecting an incident wavefront propagating through a turbulent 
medium. A particularly interesting application of this technique is in the field of astronomical 
ground-based imaging. The idea behind AO systems is to generate a corrected wavefront as 
close as possible to the genuine incident plane wavefront thanks to a deformable mirror (DM). 
An AO system is also composed of a wavefront sensor measuring the resulting distortion of 
the collected wavefront after correction by the DM. Based on these measured signals, the 
voltage applied to the piezoelectric actuators is computed in order to reshape the mirror. 
The tilts (first order modes) of the wavefront are corrected by a first mirror. Then, the 
DM is part of the control- loop for the correction of higher-order modes of the wave front. 
Different types of sensors (curvature sensor, pyramid wavefront sensor) may be used to esti- 
mate the distortions affecting the incoming wave-front but the most common encountered in 
existing applications is the Shack-Hartmann (SH) sensor. There also exists different type of 
deformable mirrors and we choose to study the case of the most common one. For additional 
details on basic principles of adaptive optics, see [1], e.g.. 

This paper is devoted to the design of control laws for an adaptive optics system formed 
by a bimorph mirror and a Shack-Hartmann sensor (see Figure 1). Most often, the existing 
adaptive optics systems use static models and very basic control algorithms based on frequent 
measurements of the influence of each actuator of the mirror to each output of the SH. 
This allows the computation of an interaction matrix gathering the corresponding influence 
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Figure 1: An adaptive optics system - The control loop consists in a SH sensor analyzing the incoming 
wavefront and a layer of piezoelectric sensors giving the precise position of the bimorph mirror, both of them 
allowing the calculation of the appropriate command of the DM in order to recover the genuine wavefront. 



functions. Here, our goal is to consider the design of an adaptive optics system control 
loop from a modern automatic control point of view as in [2] and [3]. This means first that 
dynamics of the different elements involved in the control-loop have to be taken into account. 
In particular, a specific dynamic model for the DM is proposed for control purpose (as already 
presented in [4] and see also [5]). Secondly, a state-space model of the turbulent phase, built 
from its frequency domain characteristics, is defined [6]. 

The main contribution concerns the infinite dimension setting introduced in this paper. 
More precisely, while in the literature, only static finite dimensional models are considered, 
a model based on a particular partial differential equation (pde) is used for the DM. We 
believe that our point of view matches well with the reduction of the size of the actuators 
and the significant augmentation of their numbers in many devices, as in AO for Very Large 
Telescopes. 

In reference [7], a thin elastic plate model of a deformable bimorph mirror is derived. 
This model is based on a periodic distribution of embedded piezoelectric patches that may 
be used as sensors or actuators. The idea is then to elaborate a robust control strategy based 
on modern control tools for distributed parameter systems [8]. Moreover, in contrast to [9] 
and [2], we do not need to compute any interaction matrix modelling the relation between 
the input on the piezoelectric patches attached to the mirror and the output given by the 
Shack-Hartmann sensor. The interaction matrix can be seen as a static model of the mirror 
whereas a more general dynamical model of the mirror is used here. 

For the sake of clarity of this study, we emphasize here the main informations about 
the frame we choose for our modeling of robust control of an AO system. We consider a 
continuous time state-space model of an AO loop (as in [5] and instead of a discrete one in 
[2] and [9]) and without delay. In practice an AO system uses discrete wavefront sensing 
data with inherent temporal delays and of course it is possible to derive a discrete time 
extension of our model but it is not our point here, even if we recognize that the perfor- 
mance will somehow be affected. Our contribution relies mainly on the new pde model of 
the DM and we aim at using the Hoo control theory for infinite dimension setting in order 
to recover at least similar performance as the one of LQG control for a standard model of 
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the DM (see [9]). One should notice that our model depends only on a few physical parame- 
ters (such as the density, the stiffness... see the Bimorph mirror model subsection below for 
more details), parameters that could be considered as uncertain quantities the control law 
should take into account. Therefore, we do not need either to compute an interaction matrix 
(which is more and more complicate to compute when the number of the sensors and of the 
actuators increases as for Very Large Telescopes), or the inverse of this interaction matrix [1]. 

The control problem is solved using an control setting. The first motivation is that 

Hoc control theory provides intrinsic properties of robustness while optimizing on the worst- 
case performance. Another motivation is the multivariable nature of the control problem 
involved in adaptive optics system design [3]. Current adaptive optics control systems use 
decoupling modal control to rewrite the original problem as several decoupled single input 
single output control problems. Because Hoo control framework may easily handle a mul- 
tivariable dynamic model of the bimorph DM in the synthesis process, the obtained robust 
controller outperforms usual static control approaches of the literature. In addition, the use 
of Hinfinity controllers induces, in general, some robustness properties of the closed-loop 
while HijLQG controllers (privileged in general, see [9]) lead to improvement of the perfor- 
mance but with no robustness guarantee (sec [10]). So far, we do not claim to have solved 
the complete problem of AOS synthesis (with delays and limitations of performance intro- 
duced by sampling) but we think that this new setting will probably address fundamental 
issues encountered in the very large telescopes context. This work is meant to illustrate the 
realizability of such an approach on realistic instances of AOS Design. 

The outline of the paper is the following. First, the adaptive optics control system is 
described (see Section 2) through the presentation of the models of the bimorph mirror and 
the turbulent phase. The third section is dedicated to the robust i?oo control setting in the 
infinite dimension framework and its formulation in our particular case. The last section 
contains the description of the truncated model and the numerical results. 



2 The adaptive optics model 

The bimorph mirror is composed of a purely elastic and reflective plate equipped with piezo- 
electric actuators (in order to deform the shape of the mirror) and piezoelectric sensors (to 

measure the effective deformation). A Shack-Hartmann sensor then analyzes the resulting 
phase 0res of the wavefront, after reflection in the deformable mirror of the turbulent phase 

?!'tur- 

Different types of disturbances have to be faced with: Wmod represents unstructured uncer- 
tainty (neglected dynamics) affecting the model, Wpiezo and w^^ are noise signals respectively 
attached to piezoelectric and Shack-Hartmann sensors. Finally, (jf)tur is the turbulent phase 
of the wavefront introduced by the atmospheric perturbation. 

We denote by e = e{r,9,t) the transverse displacement of the circular mirror at point 
of polar coordinates (r, 0) and time t, while A is the light wavelength. The corrected phase 
produced by e is then given by ^cor = leading to a resulting phase: 

47r 

<Pres = + ^tur (1) 

The optic sensor's output, computed by Shack-Hartmann sensor is: 

47r 

VsH = — + <?^t"'- + "^^SH (2) 
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where c is a modelling parameter of the perturbation. 



Finally, wc note that the control input is the voltage u applied to the piezoelectric ac- 
tuators and the corresponding piezoelectric output is the voltage ype measured with the 
piezoelectric inclusions used as sensors (see equations (3) and (4) below). Indeed, in com- 
parison with many other devices, where the only information used to compute the voltage u 
comes from the wavefront analyzer, the additional possibility of measuring the deflection of 
the mirror through a layer of piezoelectric sensors (see Figure 1) is considered here. 

It is recalled that the goal of the adaptive optics control system is to minimize the resulting 
phase of the wavefront using Shack-Hartmann measurements. 

Bimorph mirror model 

To obtain the model of a bimorph mirror (see an outline in [4]), we consider three different 
layers. One is purely elastic and reflective, the second one is equipped with piezoelectric 
inclusions used as actuators, the third one is equipped with piezoelectric inclusions used as 
sensors. The heterogeneities are periodically distributed. In reference [7], the authors derive 
the following dynamical model of the mirror (a partial differential equation with respect to 

p dtte + QiA^e + Q2e = dsiAu + pbwmod (3) 

with the initial conditions e(r, 6,t = 0) = eo(r, 0) and dte{r, 6,t = 0) = ei(r, 6). The voltage 
ype computed by the piezoelectric sensors is given by 

Vpe = 631 Ae + dWpe- (4) 

The following notations are defined: 

• (r, 9) are the spatial coordinates of a point of the disk fl of radius a and t is the time; 

• A is the Laplacian operator and for a general function v{r, 9) in polar coordinates 

^ dv 1 dv 1 dv 

• M is the voltage applied to the inclusions of the actuator layer; 

• p is the surface density, v is the Poisson ratio of the mirror's material, Qi is the stiffness 
coefficient and Qi is a correction coefficient; 

• 631 and £^31 are proportional to the piezoelectric tensor coefficient (isi (for more physical 
details see [11]); 

• h and d are linear applications on appropriate spaces; 

• w'mod and are unknown perturbations modelling the model errors of the plate 
equation and the measurement noise of the piezoelectric output. 

The boundary conditions are those of the free edges case (VLT and the experimental 
device SESAME, see Subsection 4.2): 
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Turbulent phase model 



In order to complete our optics system model, we need to develop a model of the turbulence 
phase. 

A usual representation of atmospheric phase distortion is made through the orthogonal 
basis of Zernike polynomials because the first Zernike modes correspond to the main optical 
aberrations. An infinite number of Zernike functions is required to characterize the wavefront, 
but a tnmcatcd basis is used in general for implementation purpose. Note that a 14-th order 
approximation contains 92% of the phase information, without taking into account the piston 
mode which represents the average phase distortion [9]. The tip/tilt modes are not part of 
our modelling of the turbulent phase because of their correction by a dedicated mirror. Wc 
will therefore work with the 12 first modes of Zernike given in reference [12] and recalled here 
(see Table 1), excluding the three first ones. 
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Table 1: First 15 Zernike Functions 



The turbulent phase (?!)tur is approximated as follows: 



Nz 



^tur(r,0,i)«='l]'/'i(i)^i(r-,( 



i=4 



where Nz > 15. is the i-th Zernike function and for all i, (piit) is a random time- varying 
coefficient corresponding to the projection of ^tur on Zi. 

= <tl 



H{3w) 



Figure 2: Shaping filter generating (j) - The turbulent phase </> is modeled through a hnear shaping 
filter of transfer function H from the noise w 

To build a state-space representation of the turbulent phase, i^itur is modelled as the 
output of a linear shaping filter (illustrated by Figure 2) of the form : 



F^ + Gw 



(6) 
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where = {4>4, ■ ■ ■ , 4'Nz)y ^ = (^''4) ' ' ' iWNz)i P ^^'^ G are two time- invariant square ma- 
trices of {Nz — 3)-dimension and w is a stationary zero-mean white gaussian noise. (f)tnv is 
therefore a stationary process. 

In order to compute F and G, the results presented in [6] and based on the Kolmogorov 
theory of turbulence and associated approximations in the frequency domain are used here. 
They confirm similar results proposed in [13] and complete the stiidy of frequency domain 
behavior for each Zernike coefficient. Each Zernike function's spectrum are characterized by 
a cut-off frequency whose heuristic expression is given by: 

/,,~0.3(n, + l)^ (7) 



where rij is the radial order of the Zernike number i 
diameter of the circular aperture of the telescope. 

The random process is supposed to be composed of Nz 
processes. For i = 4 • • • Nz, we have: 

Up) 1 



V is the average wind-speed and D the 
3 decoupled first-order Markov 
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Table 2: Atmospheric phase distortion state-space model with the average wind-speed V = 
-1 D 

9m s , — =8 and the wavelength A = 550nm 



The matrix G is obtained from the steady-state Lyapunov equation verified by the cor- 
relation matrix P^(oo): 

GG' = -{FP^{oo) + P4,{oo)F') (9) 
A closed- form expression for the spatial covariance matrix is given in [12]. 
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where T is the Gamma function and ro is the Fried parameter (corresponding to the strength 
of the turbulence [1]). Table 2 shows the non zero entries of the matrices F and G for 

V = 9m s~^ and — = 8 (as in reference [2]). 



3 Robust Control Results 

The point of this section is to prove that the new model we propose for AO systems is valid 
for an iJoo-control study. One of the difficulties comes from the infinite dimensional setting. 
For a survey of the ifoo-control theory for the infinite-dimensional case, the interested reader 
may have a look at [14] or [15] for the state-feedback case and [8] for the output-feedback case. 
The main results arc a generalization of finite-dimensional regular _ffoo-control problems (see 
for instance [10]). In particular, the solution will be given in terms of the solvability of two 
coupled Riccati equations. 

The linear infinite-dimensional model derived from the partial differential equations pre- 
sented in Section 2 has to fit in the following standard formalism of measurement-feedback 
control 

x' = Ax + Biw + B2U 

z = Cix + D12U (P) 
y = C2X + D21W 

where x is the state of the system, u is the control input, w is the disturbance input, y is the 
measured output and z is the controlled output. 




K 



y 



Figure 3: Closed-loop system 
Therefore, we introduce the following notations: 

• the state vector x = (e, die, (piui)'^ where e is the transverse displacement of the plate 
and (ptur is the projection of the turbulent phase on the first Zernikc modes; 

• the exogenous disturbance inputs vector w = (wmod,w'sH)'"^tur,w'pe)^ gathers the dif- 
ferent perturbation signals (uncertainty affecting dynamics of the model and of the 
turbulence phase, noise vectors of the wavefront analyzer and of piezoelectric sensors); 

• the control inputs vector u is the voltage applied to piezoelectric patches; 

• the measurement outputs vector y = (ype,2/sH) composed with the piezoelectric and 
the wavefront analyzer measured outputs; 

• the controlled outputs vector z = (cprcs, u) contains an optical part (the resulting phase, 
see (1)) and the control input vector u. 
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The aim is to find a dynamic measurement-feedback controller K ensuring that the in- 
fluence of w on 2; is smaller than some specific bound. The corresponding standard block 
diagram is given by Figure 3. 

The controller is assume to have the following form: 



u 



Mp + Ny 
Lp + Ry 



(K) 



where M is the infinitesimal generator of a Co-semigroup on a real separable Hilbert space 

and N, L and R are bounded linear operators. With this controller, the closed-loop system 
can easily be derived and defines a bounded linear map Sk such that z{t) = {SKw){t). Its 
bound is denoted ll^i^rHoo- 

The control loop defining the adaptive optics system is sketched in Figure 4. If we gather 
the different equations describing the system, namely (1), (2), (3), (4) and the forthcoming 
equation (11) (corresponding to (6)), we get 



due + QiA^e -I- = d^iAu + bpw^od 

dt<l>tui = -^^tur + GWtui 
0res = 0tur -e 



Vpe = 631 Ae -r u,uvpe 
ysH = -^tij 



A 

47r 

— e + cWg 



(10) 



Actually, in order to have an unified infinite dimensional modelling of the adaptive optic 
system's state, we described the model of 0tur from equation (6) as follows: 

• "^^tur and Wtur are the reconstruction of (f) and w on Nz — 3 of the first Zernike modes, 
such that 

Nz Nz 
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what leads to the turbulent phase model given in (10) 



't<^tur 



Gwtu 



(11) 



where is the Hilbert space of square integrable functions and C{X) stands for the set of 
linear applications on X. 

Thus, the operators defining the standard form P are built from (10) 
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Figure 4: Standard model for adaptive optics system control loop 

The appropriate functional spaces associated to the infinite-dimensional model are now 
precisely defined. With the boundary condition (5), we consider the state space (the mirror 
J7 is a disk of radius a) 

X = Hl^{n) X L^{Vl) X L^{Vl) 

= {e e H^{fl),e satisfying (5)} x {L^in)Y 

the input spaces U — H'^{fl) n i?o(f^) and W = (i^(ri))^ and the output spaces Y = Z = 
(L^(ri))^, where Hq and are the Sobolev spaces 

H^{n) = {/ e L^n)/ Vz = 1,2, dj e L^{n), f\on = o} 
H^(n) = {fe L^m/ vz,j - 1,2, dj, d.djf e L^n)} 

This model satisfies all the assumptions of the main theorem of reference [8]. We give 
here a simplified version of this result: 

Theorem 1 [8] Let 7 > 0. There exists an exponentially stabilizing dynamic output-feedback 
controller of the form (K) with ||S'a'||oo < 7 if and only if there exist two nonnegative definite 
operators P, Q Cz C{X) satisfying the three conditions 

(i) y-x £ D{A), Px e D{A*), 

{A*P + PA + P{-f-^BiBl - B2B;)P + ClCi)x = 

and A + {'j^'^BiBl — B2B2)P generates an exponentially stable semigroup, 
(a) Vx e D{A*), Px e D{A), 

{AQ + QA* + Q{j-^CIC\ - C^C2)Q + B,Bl)x = 
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and A* + (7 ^Cj* Ci — C2C2)Q generates an exponentially stable semigroup, 
{Hi) 

r„{PQ) < 7^ 

where ra-{PQ) stands for the spectral radius of PQ. 
In this case, the controller K given by {K) and 

M = A+ {j-^BiBl - B2B*)P 

-Q{I-j-^PQ)-^CIC2 
N = -Q{I-^-^PQ)-^C^ (12) 
L = B^P 
R = 

is exponentially stabilizing and guarantees that we have \\Sk\\oo < 1, *e 

||<;^res|U2(a) + ||u||l,2(o) < 7lKII(Z,2(Q))*- 

Finally, if the solutions to the Riccati equations exists, then they are unique. 

Upon additional assumptions that are not detailed here, the main point is to prove that 
A is the infinitesimal generator of a Co-semigroup on the real separable Hilbert space X. 
Actually, if we consider the unbounded linear operator 



Ai : V{Ai) 





where 

V{Ai) = {eo e H^{n),eo st (5)} x H'^{n) x L'^{n), 
then one can prove that Ai is dissipative on X. Indeed, we prove that for all x G X, 

{A-ix,x)x < 

using the following scalar product on H^^{Q) in cartesian coordinates {xi,X2) G ^, as sug- 
gested in [16]: 



<u,v >H^^in) 
= / AuAv-{l-u)(—2^ + 



2(1 ( ^"^^ \ dQ 

\dx1dx2 dxidx2 J 

Moreover, one can easily check that Ai is also self-adjoint and onto. Therefore, from Lumer- 
Phillips' Theorem (sec [17], p. 15), Ai generates a continuous semigroup of linear contractions 
acting on X. And finally, since A is the sum of Ai and of a linear operator bounded on X 
(as is assumed to be bounded, like F), the proof is complete (see [18], p. 40). 

Of course, from a numerical point of view, we need to get an appropriate finite dimensional 
model. 
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4 A truncated model for numerical design 



4.1 Truncation 

The corresponding finite dimensional model can be presented as : 

x'jq = Anxn + BinWn + B2NUN 

zn = Cinxn + D12NUN (13) 

UN = C2NXN + D21NWN 

where the operators of system (P) have been replaced by real- valued matrices computed on 
truncated hermitian basis. We denote by Nb the number of eigenfunctions of operator we 
consider and by Nz the number of Zernike modes used to describe (ptui- Then, xn G M.'^^b+Nz 
is the state vector, Wn & ]^'^Nb+2Nz jg ^j^g exogenous perturbation vector, un G is the 
control vector, zn S R^^+Nz jg ^j^g controlled output vector and j/jv S R^^+Nz jg ^jjg 
measured output vector. The matrices An, Bin, B^n, Cin, D12N, C2N and D21N are of 
appropriate dimensions. 

In order to compute these objects, we still consider the case of a circular bimorph mirror 
which is free at all the boundary (this is also the case of the mirror considered in Section 4.2 
below). The eigenvectors of operator 

P P 

are given by, for all {k,j) G N^, 

Lkj{r,6) = akj (^Jk (-^^ + Ckjh (^"^^^ cos(fc6i) 
Mkj{r,6) = akj [jk + Ckjh (^)) ^H^O) 

where {r,9) are the polar coordinates of a; e fi, Jfe and Ik are, respectively, ordinary and 

j — ^ the corresponding 

eigenvalues. The family 

{Lfe,,Mfc,,(fc,j) eN^} 

is an Hilbertian basis of H^^{Q). The dimensionless coefficients Xkj and Ckj depend on the 
boundary conditions while Ofcjis computed using a normalization condition on the eigenvec- 
tors (see [19] for further details). In what follows, we consider the case of Poisson ratio 
v = 0.2 corresponding to the material the mirror is made of. Once a maximal azimuthal 
order is given (here fcmax = 5) the modes are classified according to increasing Xkj and one 
has the values gathered in Table 3. 

The sequence of functions Lkj and Mkj need to be re-ordered. They are now denoted by 
Bn and follow the increasing values of Xkj, alternating cosine and sine and eliminating the 
null eigenvectors Mqj. Therefore, 

yxeX, x= otiBi{r,0) 

neN, i>l 

where (a„)„>i is a sequence of real numbers satisfying X^^g^j „>i < 00. 

In reference [20], one can find that this basis (i?„)„eN with free boundary conditions is not 
orthogonal in L^(f2). However, numerically, we can prove that this basis is nearly orthogonal, 
indeed lots of scalar products in L'^{Q.) are null and the others are small (10~^) in comparison 
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Table 3: Coefficients of the eigenvectors L^j and Mkj- 



with unity. So, for more numerical faciUties, we will use the scalar product in Li^ip) rather 
than in Hl.{9.). 

Given Ng and Nz G N, wc compute Ajy, Biat, B2N1 CiN, C2N1 D12N and D21N using 
the "Bessel" truncated basis {Bq, Bi,. . . , B^^} and the Zernike one {Zq, Zi, . . . , Zj^^}. 

We make analogous assumptions for the tuning parameters 6, c and rf, i.e. h = diagj(6j), 
c = diagj(Q) and d = diagj(di) where (&i)ieN, i>i, {c-i.)iefL t>i and (di)^^^, i>i are sequences of 
real numbers. We recall that these coefficients are weighting functions defining the respective 
weights of the disturbance signals and the choice of diagonal matrices corresponds to an 
assumption of decoupling between the different modes. 

Futhermore (pres is expressed on Besssel functions, so we need to estimate a projection 
matrix to define 0tur with Bessel spatial coordinates. We note Q this projection Nb x Nz- 
dimension matrix. Thus, the computed equation becomes: 

4n 

4'res,i = + QijBj+2 

We denote by each null matrix with the appropriate dimensions so that each following 
matrix makes sense. We get 
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where cof = — f — | + — and (•, •) is the usual scalar product in L'^(fl). 







4.2 Numerical results 

In this subsection, numerical simulations are proposed. To get more realistic results, the 
experimental device of the project SESAME of the Observatoire de Paris is considered. This 
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VV IIIQ spCCQ 


V — y Ills 


Diameter of the pupil 


D = 10"^ m 


radius of the mirror 


a = 25 X 10"^ m 


mirror's stiffness coefficients 


Qi = 84 Nm, Q2 = 11.25 x 10* Nm"^ 


mirror's surfacic density 


p ~ 16.3 kg.m^'' 


piezoelectric coefficients 


dsi = -0.0044 NV-i, 631 = -5.60 x 10^ Vm 


wave length 


A ~ 550 nm 



Table 4: Physical parameters for the numerical simulations 



experimentation uses a biniorph mirror with a distribution of 31 piezoelectric actuators. The 
piezoelectric inclusions are PZT patches. We use the physical constants of Table 4 

We simulate only the 12 modes which follow the tip/tilt. The performance of the control 
system is evaluated by considering the spatial norm ||.||i2 of 4>ros compared to ||</>tur||L2(-Q'): 

For identical random initial conditions and taking the respective weights of the disturbance 
signals such that bi = 0.001, Ci = 0.002 and di — 0.003 for all i, we obtain the results 
represented in Figure 5. 

0.045 



0.04 
0.035 
I 0.03 
0.025 

0.02 
0.015 

0.05 0.1 0.15 0.2 0.25 

time (s) 

Figure 5: Time-evolution of |l0tur|lL2(a) (solid line) and ||?!'ros||L2(n) (dashed line) 

Using Monte Carlo simulations, the ratio between temporal average of H^turllL^j-Q-) and 
||0rcs||L2(Q) is near to 1.91 which represents a phase distortion attenuation of the reflected 
wavefront of 48%. In addition one should recall that this result does not take into account 
the tip/tilt correction. Even if these results are of the same order of magnitude as those 
presented in [9], which cannot be considered as completely satisfactory when considering 
usual results on real experiments, they clearly demonstrate the feasibility of the proposed 
approach. The possible degradation of such a performance induced by the delay in the loop 
and the discretization of the control law for implementation purpose could darken the picture. 
It must be recalled that this apparent loss of performance is mainly due to the tuning of the 
trade-off between robustness and performance that is inherently encountered in closed- loop 
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feedback design. Numerous improvements have still to be considered as presented in the next 
conclusion. 

5 Conclusion 

In this paper, a new framework to deal with the problem of adaptive optics is proposed. It 

is mainly based on an infinite-dimensional model of the dcformablc mirror associated with 
the definition of a standard model on which robust control techniques may be applied. The 
preliminary numerical experiments show a performance level comparable with the results 
of reference [9]. The main advantage of the approach suggested in this paper is that no 
interaction matrix is required to control the system. We do not pretend to outperform 
already existing AO systems but rather to pave the way for future major improvements in 
terms of robustness and efficiency of the proposed control strategics. The authors arc planing 
to take into account a model for the Shack-Hartmann wavefront sensor including a time delay 
associated with processing measurements. This will be covered in a next study. 
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